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Abstract. 

We develop further the integration procedure in the generahsed invariant for- 
malism, and demonstrate its efficiency by obtaining a class of Petrov type N 
pure radiation metrics without any explicit integration, and with comparatively 
little detailed calculations. The method is similar to the one exploited by Edgar 
and Vickers when deriving the general conformally flat pure radiation metric. 
A major addition to the technique is the introduction of non-intrinsic elements 
in generalised invariant formalism, which can be exploited to keep calculations 
manageable. 

PACS numbers: 0420, 1127 



1 Introduction 

The NP formalism [17] has been a very useful tool for obtaining and analysing 
solutions of Einstein's equations. The tetrad freedom inherent in this formalism 
has been both a boon and a hindrance; often a skillful tetrad choice gives an 
immediate simplification which sometimes leads directly to a satisfactory con- 
clusion, but at other times the original simplification proves to be cosmetic and 
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an impass is reached in the calculations. There is usually an obvious choice 
for the direction of the first null tetrad vector; sometimes the direction of the 
second null tetrad vector suggests itself; but there is usually no geometric mo- 
tivation to fix the spin and boost freedom. An additional disadvantage is that 
the exploitation of the tetrad freedom, especially that associated with null ro- 
tation and spin and boost after the first null direction is fixed, can require long 
painstaking calculations, keeping track as intertwined tetrad and coordinate 
freedoms are gradually used up in successive steps; the possibility of errors is 
always present. 

The GHP formalism [8], which generalises the NP formalism, has the spin and 
boost freedom inbuilt (as spin and boost 'weight') which means that we do not 
need to make an explicit choice of spin and boost gauge; so this is the ideal 
formalism to use when investigating spaces which pick out two null directions 
in a geometric manner. Moreover, by rescaling all quantities in the GHP for- 
malism to zero-weighted quantities, together with only one complex non-zero 
weighted quantity, we are able to seperate out the effect of the spin and boost 
gauge and essentially ignore it; working exclusively with zero-weighted quanti- 
ties gives natural coordinates and results which are more accessible to physical 
interpretation. In addition the messy calculations often associated with fixing 
the spin and boost gauge do not occur. 

However, investigations of those spacetimes where only the direction of one null 

vector is picked out geometrically, often encounter the complication associated 
with the freedom of a null rotation. This freedom can cause analogous compu- 
tational complications in a GHP analysis as the spin and boost freedom case in 
a NP analysis. 

The generalised invariant formalism (GIF) [16], [15] generalises the GHP for- 
malism by essentially building the null rotation freedom into the GIF, which 
means that the formalism is built around one dyad spinor O^, and we do not 
need to make an explicit choice for the null rotation freedom of the second dyad 
spinor la- But the price we pay is that we have to move from a scalar to a 
spinor formalism, which at first sight appears very complicated, and even un- 
wieldy. However, manipulating within the GIF is in practice (although not in 
appearance) not much more complicated than the scalar manipulations of the 
GHP formalism. 

Some years ago Held [9], [10] outlined an idealised integration method in an 
'optimal situation' using the GHP formalism. He advocated manipulating the 
GHP equations until they were reduced to a complete and involutive set of ta- 
bles involving first derivative GHP operators; from these tables Held argued it 
should be possible to read off the metric without having to perform any explicit 
integrations. Subsequently this idea has been developed and applied in a num- 
ber of ways, and it is now understood that the 'optimal situation' is when the 
manipulation succeeds in generating a table for each of four real zero-weighted 
scalars (which will become the coordinates) and a table for one complex non- 
zero-weighted scalar (which describes the spin and boost gauge) [4]. An impor- 
tant technique in this method is the repeated application of the commutator 
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equations; however, in addition, it is crucial to recognise the theoretical impor- 
tance of the commutator equations where important information resides (as well 
as in the GHP Ricci and Bianchi equations) and in order that all this informa- 
tion has been extracted it is essential that the commutators should have been 
applied explicitly to these five scalars [4]. These five tables will not in general 
be completely involutive and we would expect some additional scalars and their 
(possibly, partial) tables to ensure a complete and involutive system. 
Following the pioneering work of Held [9], [10], the most striking illustrations 
of the advantages and efficiency of this GHP formalism have been illustrated 
in a re-investigation of Petrov D vacuum metrics [1], [2] where key constraint 
equations were obtained with a fraction of the calculations required in the NP 
formalism, and in the re-investigation of the conformally flat pure radiation 
metrics [5] where new solutions were obtained, which had been overlooked in 
the complexities of an earlier NP investigation [18]. 

In a recent paper by Edgar and Vickers [7] this GHP integration method was 

generalised to the GIF [15]. Again, the method consists of manipulating all the 
equations of the formalism to construct a complete and involutive set of tables 
involving first derivative GIF spinor operators. The 'optimal situation' to be 
sought is to generate a table for each of four real zero-weighted scalars (which 
will become the coordinates), a table for one complex (non-trivially-) weighted 
scalar (which describes the spin and boost gauge) and a table for a second spinor 
1a', such a spinor (which is not parallel to the first dyad spinor Oa) should emerge 
naturally from the calculations, and can then be identified as the second dyad 
spinor La- Again, an important element in this method is to recognise that 
much information resides in the GIF commutator equations (as well as in the 
GIF Ricci and Bianchi equations) and in order that all this information has 
been extracted it is essential that the commutators should have been applied 
explicitly to these five scalars as well as to the new spinor . Once these tables 
have been found, and the new spinor identified with the second dyad spinor La, 
the problem can be reduced to a purely scalar one in the GHP formalism. 
The advantage of the GIF was illustrated by a further investigation of the con- 
formally flat pure radiation metrics [7]. In the GHP investigation a degree of 
guesswork and luck was involved in obtaining that particular choice of null ro- 
tation which fixed the direction of the second null vector in such a manner that 
subsequent calculations became manageable. In the GIF no explicit choice was 
needed, since the formalism deals only with quantities which are invariant under 
such null rotations. Moreover, the resulting form for the class of metrics was 
such that Killing vector properties were obvious and the invariant classification 
procedure could easily be deduced; such conclusions could not be made from the 
form of the metric obtained by the GHP analysis. The two metric forms differ 
only slightly in appearance, but fundamentally different in the interpretation of 
some of the coordinates. 

An intriguing aspect of this GHP and GIF integration method is that it is 
best suited to spaces without any Killing vectors. The principle of the method 
is to try to generate the four coordinates directly from the intrinsic elements 
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of the formalism, and in spaces with no KiUing vectors the coordinates will be 
generated in this manner. When Killing vectors are present there are not enough 
intrinsic elements of the formalism to generate four intrinsic coordinates, and 

so new non-intrinsic elements have to be introduced. 

The conformally flat pure radiation metrics [5], [7] provided an 'optimal situa- 
tion' and gave a comparatively simple demonstration because we were able to 
generate all four cordinates intrinsically from the elements of the formalism in 
the generic case [7]. However, there was a subclass of these spaces for which we 
could only generate three coordinates intrinsically (corresponding to the exis- 
tence of a Killing vector). We were able to identify these spaces by modifying 
the technique and using the structure of the commutators to motivate a coor- 
dinate choice from outside the intrinsic elements of the GIF; we then explicitly 
confirmed that this choice was compatible with all the other equations. 
This technique of introducing a new non-intrinsic element into the computations 
— motivated by the property that it automatically satisfies the commutator 
equations and implies no new constraints — is essential when dealing with 
spaces in which Killing vectors are present; however, this technique can also 
be used to try and obtain simpler tables when the direct method generates 
complicated ones. 

So far, the conformally flat pure radiation metrics are the only class of metrics 
which have been determined explicitly by this GIF procedure [7]. The GIF 
formalism was initially designed to be used to investigate the invariant classifi- 
cation of spacetimes; in particular of those metrics where the second null tetrad 
vector has no simple geometric links, e.g. Petrov type N vacuum spaces [14]. 
However, it is clear that the GIF can also be used for obtaining solutions of 
Einstein's equations for metrics with these properties, but before attacking very 
difficult problems like the whole class of vacuum Type N twisting metrics, it is 
necessary to build up more experience of using the GIF. 

Therefore, in an effort to better understand how the GIF integration procedure 
works in practice we shall investigate a class of Petrov type N pure radiation 
metrics. These metrics are close enough to their conformally flat counterparts 
in [7] for us to get some hints of how to tackle them in the GIF, but sufficiently 
different to present some new challenges; in particular we expect a richer Killing 
vector structure, and hence a less direct method. Unlike in [7] where we followed 
a direct method, and only brought in one coordinate candidate from outside the 
formalism at the last step, in this approach we will introduce a new element from 
the beginning. These metrics are well known, especially in their familiar Kundt 
form [11]. So the purpose of this paper is not to find new solutions, but rather to 
rederive a known class by developing new techniques within the GIF procedure; 
in particular to demonstrate how introducing non-intrinsic elements can keep 
calculations simple and manageable. We are further motivated to demonstrate 
the power of the GIF approach for these spaces because attempts to make an 
investigation of these spaces simply and efficiently in the GHP formalism have 
been imsuccessful due to the problem of the null rotational freedom of the second 
tetrad vector; instead a very detailed and involved calculation has been needed 
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in NP formalism and in a hybid approach combining GHP with NP ideas [6]. 

In the next section wc summarise the GIF [15], with special reference to dif- 
ferential operators. In Section 3 we specialise the GIF to the class of spaces 
under discussion, and in Section 4 we carry out the GIF integration procedure 
for these spaces. Section 5 gives a summary. 

2 The Differential Operators 

In the GIF the role of the spin coefficients k, a, p and r is taken up by spinor 
quantities K, S, R and T given by 

K = K 

Sa' = ooa' - kZa' 

Ra = pOA - KtA 
TaA' = TO AO A' - POaIa' - CTIaOa' + klaIa' (1) 

Under a transformation of the spin frame given by 

o^^Ao^ L'^^X-h^ + ao^ (2) 

these transform as 

K A^AK 

Sa' ^ A^S^, 
Ra ^ X^XRa 
Taa' ^ X^Taa' (3) 

They are therefore invariant under null rotations and have weight {p, q} under 
spin and boost transformations given by 



K : 


{3,1} 




S : 


{3,0} 




R : 


{2,1} 




T : 


{2,0} 


(4) 



The role of the differential operators P, d, P' and d' is played by new differential 
operators P, d, I>' and d' which act on properly weighted symmetric spinors 
to produce symmetric spinors of different valence and weight. These operators 
may all be defined in terms of an auxiliary difi'erential operator Vaba' b' which 
is defined by 

'DabA'B'VCi ■..CnC[...C'^, 

= OaOa'"^ BB'11Ci...CnC[...C'^, 

-(POa'Vbb'Oa +^OA'^BB>OA')r)Ci...CNC[...C'^, (5) 
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(^)aCi...CnA'C[. 




= 2?^BA'B'J?Ci...CivC;...i 






sym 


{9ri)ACi...CNA'B'C[. 




= '^o'^VabA'B'Vc^...CnC[...C'^, 






sym 


{0'rj)ABCi...CNA'C[. 




= 1^ABA'B'r)Ci...CNC[...C'^, 






sym 


{^'v)aBCi...CnA'B'C[. 




= ^T^ABA'B''nCr...CNC[...C'^, 



where r) has weight {p, q}. 

We will need to know the result of contracting P'jj with o and 6. We can write 
equation (5) in the form 

1^aba'B'Vcu..CnC[...c'^, 

= (i''VCi...CNC{...C'^,)^AOBOA'OB' 
-(^'^Ci...C«C;...C;^,)OAOBOA'ti3' - id'nCi...CNC[...C'^,)(^Al'BOA'OB' 
-i^nCi...CNC[...C'^,)OA''BOA''iB' 
+ {PLaOa'TbB' +<iOAlB''TB' 3)1101. ..Ci,'C[...C'^„ (6) 

where P', d', d and P are the ordinary GHP operators applied to spinors. The 
new operators are obtained by contraction with o and 6, and symmetrizing. 



(8) 



(10) 

where indicates symmetrization over all free primed and unprimed indices. 

sym 

In the case of a scalar field this gives 

{^V)aBA'B' = (P'?7)OaObOa'Ob' - {d'r] - qfri)OAOBO(A~l'B') 

-{drj - pTri)0(A''B)OA'OB' + i^V - PPV - QP-ii)0(aI-b)0(a'Ib') 

+ {j)Kt.Al'BO(A'lB') + qKO(At'B)iA'lB' 

-patAt-BOA'OB' - qa-OAOBl'A'tB')'n (11) 

Contracting (11) with 6^' gives 

{V'r])ABA'B'0^ = i{{d'r])ABA' - qifOAOBOA' - pO{Al'B)OA' 
-UOaObIa' + KO(AiB)iA')^} 

= k{{d'r})ABA'-qTA'{AOB)ri) (12) 
or in the compacted notation 

(P'r,)-6=i{(5''?)-aTr?} (13) 
Similar calculations give 

(P'r/)-o=i{(«?y)-pTr?} (14) 
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{d'v)-o=^{{Prj)-pKrj} 



(15) 



(dv) ■ 6 = i{{p7]) - qKrj} 



(16) 



(P'ry) • • O = iiiVr]) - pRr] - qR?]} 



(17) 



For a spinor f] the above contractions become more complicated. For example 
for a valence (l,0)-spinor jj^ of weight {p, q} we get 



Although the definition of the differential operators is quite complicated, the 
fact that they take symmetric spinors to symmetric spinors means that one can 
write down the equations in an index free notation. 

The Ricci equations, Bianchi equations and the commutators for the general case 
are given in [2] . This complete system of equations is completely equivalent to 
Einstein's equations, and to find solutions to Einstein's equations this system 
will therefore have to be completely integrated. However, in view of the more 
complicated nature of the operators in this formalism, some of the information 
which resided in the Ricci equations in NP and/or GHP formalisms is contained 
within the commutators in this formalism; in particular these commutators 
contain inhomogeneous terms explicitly dependent on the weight and valence of 
the spinor on which they act. To extract all the information in the commutators 
we need to apply them to, [15] 

(i) four functionally independent {0, 0} weighted real scalars, 

(ii) one {p, q} weighted complex scalar where p ^ ±q, 

(iii) one valence (1,0) spinor, of weight {—1,0}. 

Of course, we can extract all the information by applying the commutators to 
different (but essentially equivalent) combinations of these scalars and spinor; 
however the particular choices above are best suited to our integration proce- 
dure since the four {0,0} weighted real scalars will become the coordinates, 
the complex scalar is given by the gauge field PQ, while the spinor will be 
identified with the second dyad spinor t^. 

3 The equations 

We restrict attention to the Petrov type N pure radiation spaces within the 
Kundt class of spacctimes (with a non-expanding and non-twisting null congru- 
ence); this means that when we choose Oa to be aligned with the propagation 
direction of the radiation that the only remaining parts of the Riemann tensor 
are $22 and ^'4 which we simplify to $ and ^ respectively, remembering that 
$ is real; also p = = a = k. In GIF, the Ricci spinor takes the form 



(P'jj) • o = i{V'{r, ■ o) + {d'r,} - (P - l)Tn} 



(18) 



^ABA'B' = $OaObOa'Ob' 



(19) 
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where $ is a real scalar field of weight {—2,-2}. and the Weyl spinor takes the 
form 



*ABCD = *OaObOcOd (20) 

where \1/ is a complex scalar field of weight {—4, 0}. The spin coefficients in GIF 
satisfy 

K =0 
S =0 

R = (21) 

but 

Taa' = tOaOa' (22) 

Notice that t, 'I'4 and $22' arc all invariant under the group of null rotations 
so that they can be used instead of their GIF spinor equivalents; this gives a 
considerable simplification in the GIF notation. 

The GIF equations are: 

(i) GIF Ricci equations: 



Pr = (23) 

dr = (24) 

d'r = TT (25) 

(ii) GIF Bianchi identities: 

M = (26) 

P$ = (27) 

d^-d'^ = r*-r$ (28) 

(iii) GIF commutators (applied to an invariant spinor rj): 

(W -irp)r, = {Ta + T8')rt (29) 

{Vd dP)r] = (30) 

idd' d'd)r] = (31) 

(P'd-dP')T7 = -rP'jj - $(J7 • o) - ^(j? • o) (32) 



These GIF equations contain all the information for the Type N pure radiation 
metrics. We emphasise that we assume throughout that r ^ 0. 
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4 The integration procedure 



4.1 Preliminciry reEirrangement. 

The spin coefRcient r will supply one real zero- weighted scalar (rt) and one real 

{1, — l}-wcightcd scalar (t/t). However to keep the presentation of subsequent 
calculations as simple as possible, it will be convenient to rearrange slightly, and 
following [7], use instead the zero- weighted real scalar, 

A=^=. (33) 

and the complex scalar P, 

(34) 

where P is a scalar of weight {1,-1}, and PP = |. 

We are assuming r 7^ 0, and so A and P will always be defined and different 

from zero. 

These choices give the two simple partial tables 

PP = 
dP = 

a'P = (35) 

and 

M = 

dA = -p 

d'A ^ -P (36) 

Wo could now complete these tables with the fourth operator P' by introducing 
unknown spinors and then applying each of the commutators to the completed 
tables of P and A respectively; unfortunately, the calculations soon become 
quite complicated. 

Alternatively we could look for the other part of the weighted complex scalar by 
using expressions such as the {—1, — l}-weighted scalar used in [7]. Again, 
calculations soon get quite complicated. 



4.2 Finding tables for the complex scalar PQ, and the 
spinor I and applying commutators. 

So instead we examine the commutators and check whether their comparatively 
simple structure suggests the existence of some simple table(s). In particular 
we know that we require a table for a real weighted {p = q ^ 0) scalar to 
combine with P to give a non-trivial weighted scalar of weight p 7^ ±q. From 
the structure of the commutators we are motivated to consider the simplest 
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possible partial table for a {—1, — Ij-weighted scalar Q annihilated by the first 
three operators 



PQ = 
dQ = 

d'Q = (37) 

which leads us to consider the table for PQ 

P{PQ) = 

BiPQ) = 

d'{PQ) = 

P'(Pg) = -|l (38) 

where we have completed the table with some spinor I, which is as yet undeter- 
mined. (We have introduced the weighted factor ^ in the above definition for 
I simply for convenience in later calculations.) 
It follows from (13) and (14) that 

1.5 = -^{P{PQ)) ■ 6 = -^d'iPQ) = (39) 



1.0= -q{^'{PQ)) ■ o = -Q{d{PQ) + 2t{PQ)) = -1 (40) 
Hence I is a (1,0) valence spinor, and from 

i^'iPQ))ABA'B' = -|I(AOb)6a'6b' (41) 

we conclude that its weight is {—1, 0}. 

So now we have to apply the commutators to the table for {PQ) which yields a 
partial table for the spinor I; the complete table can be written as 

W = 

01 = 

d'l = 

l/l = (42) 

where the spinor W is as yet undetermined. (The factor ^^j- is again just to 
improve efficiency of presentation.) It follows from (18) that 

W-o=-^(P'l).o = (43) 
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Hence 



where W is a (2,0) valence spinor of weight {2,0}, and is a zero- weighted 

complex scalar. 

We next have to apply the commutators to the previous table for I after the 
substitution (45); we obtain a partial table for the zero- weighted complex scalar 
W, 

= 
8W = -^—^ 

d'W = (46) 

When the commutators are applied to this partial table we obtain 

= (47) 
1)$ = (48) 

d^-d'^ = T*-T$ (49) 

which are precisely the three Bianchi equations (26)- (28). Hence our guess for 
a table for Q is completely compatible with the commutators and all the other 
equations. Furthermore, checking the compatibility of our choice of table for 
PQ led to a table for the spinor I, whose compatability we have also checked 
via the commutators. 

So we have obtained two of the core elements required in our analysis — a 
weighted scalar PQ and a new spinor I which is not parallel to O — and used 
them to extract information from the commutators. 

4.3 Finding tables for two coordinate candidates and ap- 
plying commutators to each. 

Wc can now use Q and I as defined above to complete the table for A, respecting 
spinor valences and weights and identities (13) and (14), and obtain 

PA = 

dA = -P 

d'A = -p 

V^A = PI + Pi+^N, (50) 
where A'' is a real zero-weighted scalar, which is as yet undetermined. 
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The application of the commutators to this table determines a table for the 

scalar A'', 




(51) 



where L is a real zero-weighted scalar, which is as yet undetermined. 

The application of the commutators to this table determines a partial table for 

the scalar L, 



At this stage we will review what we have obtained so far. We have tables for 

I, P, Q, A, N respectively, to all of which wc have applied the commutators; this 
has yielded the further partial tables for L and for (complex) W. Since A and 
N are real zero-weighted scalars to which we have applied commutators they 
arc two obvious candidates as coordinates. (Of course it is necessary to confirm 
that these scalars are functionally independent, a fact which would seem obvious 
from the structure of the right hand side of their respective tables; however, we 
will take care to confirm this explicitly when we translate these tables into the 
purely scalar formalism.) 

If we wish to adopt either of the three remaining real zero-weighted scalars {L 

or the real or the imaginary part of W) as coordinates, it would be necessary to 
complete the two tables and then apply the commutators to each; calculations 
then begin to get complicated. 

4.4 Finding tables for two more coordinate candidates and 
applying commutators. 

Alternatively, we can go back to the commutators to see if they suggest simple 
tables for the remaining coordinates 

We begin with the very simple table for a zero-weighted scalar T, 



PL = 
dL = PW 
d'L = PW 



(52) 



PT = 

dT = 

d'T = 













Q 

A 



V'T = 



(53) 
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and wc easily confirm that this tabic is internally consistent, as well as compat- 
ible with the commutators and with all the other equations. 
Next we try the slightly more complicated table for a real zero-weighted scalar 
B, 

PB = 

dB = -iP 

d'B = iP (54) 
V'B = i{PI -~PI) 

and again confirm that this table is internally consistent, as well as compatible 
with the commutators and with all the other equations. 

So we now have four coordinate candidates A,N,T,B, but it is necessary to 
confirm that these are functionally independent before we can adopt them for- 
mally as coordinates. Since these four scalars are zero-weighted, we can easily 
rewrite their tables using the NP operators. We note that there are also the 
tables for P, Q and I, but, since these contain only gauge information, these 
tables play no role in the construction of the metric, and we do not consider 
them further. 



4.5 The tables in terms of scalar operators. 

If we identify the spinor I with the second dyad spinor t, then the four tables 
for the zero- weighted coordinate candidates T,A,N,B can be easily translated 
into the familiar NP operators, 

DT = 
6T = 

6'T = (55) 

AT = I 

A 



DA 
5A 
5'A 

AA 



= 



-P 
-P 

In 



DN 

5N 
6'N 

AN 



1 

Q 



(56) 



(57) 
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DB = 
SB = -iP 

6'B = iP (58) 
AB = 

A simple observation of the determinant formed from the four vectors taken from 
the right hand sides of each of the four tables respectively, confirms that all four 
scalars arc functionally independent; so wc formally adopt them as coordinates. 
We note that these four tables for T, N, A, B are not strictly involutive in these 
scalars; there occur also the real scalars L and W which satisfy respectively 

DL = 
6L = PW 

S'L = PW (59) 

DW = 
SW = -=^$ 



PQ2 
PQ^ 



S'W = -=^* (60) 



4.6 Using coordinate candidates as coordinates. 

If we now make the choice of the coordinate candidates T, N, A, B as coordinates 

t = T, n = N, a = A, b = B 

the tetrad vectors can be obtained in the t, n, a, b coordinates from their respec- 
tive tables as follows: 

r = (0,-^,0,0) 

n' = ^(l,L,n,0) 
a 

m' = P(0,0,-l,-i) (61) 
to" = P(0,0,-l,i) 

and we can write out the NP derivatives as follows: 

A =«|+Li^ + (64) 

a ot on a oa 
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where it will be convenient to write ^ = a + ib. 



The metric g can now be constructed using 
to give in coordinates t, n, a, b, 



9'' = 



(65) 



/ 


-l/a 







-1/a 


L/a 


—n/a 








—n/a 


-1 





V 








-1 / 



(66) 



Prom (59) we have W{t, = L^^ where L is now a function of the three 
coordinates t, ^, ^, and so from (60) we have 



2A ■«« 



(67) 



The above analysis is dependent on the condition that none of our four coor- 
dinate candidates could be constant and still compatible with the equations; it 
is clear from considering the respective tables (55), (56), (57), (58) that this is 
the case. Therefore we can conclude that our analysis is complete in the sense 
that we have obtained the complete class of the metrics under consideration, 
and have not missed any 'special cases'. 

The metric given by (66), is in very similar form to its version in [6]. We have 
made no restrictions on the Riemann tensor components; therefore = gives 
the conformally flat special class, and = gives the vacuum special class. 
Note however that the conformally flat class is in different coordinates that in 
[7]. 



5 Conclusion 

We have described here an alternative method within the GIF for obtaining so- 
lutions of Einstein's equations and applied it to a particular class of Petrov type 
N pure radiation spaces. The basic ideas of this method were developed in pre- 
vious work by Edgar and Vickers [7] , where it was applied to obtain the general 
conformally flat pure radiation metric. That class of spaces were ideal because 
the absence of Killing vectors, in general, enabled a comparatively simple and 
direct analysis to be carried out; in the present paper such a direct analysis was 
not possible, and this is linked to a richer Killing vector structure. So, in this ap- 
plication a modified approach was needed, and we have demonstrated how the 
introduction of non-intrinsic elements, imtroduced in compatability with the 
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commutators, not only overcomes any problem created by the absence of suffi- 
cient coordinate candidates, but also can keep the calculations manageable. We 
anticipate that a major practical difficulty encountered when using this method 
resides in how to construct the second spinor I such that the subsequent tables 
will be as simple as possible; the procedure of introducing non-intrinsic elements 
gives us some control over this. 

A big advantage of this procedure resides in the fact that one has managed to 
avoid having to solve any differential equations. Another important advantage 
is that the complicated detailed calculations needed to keep track of coordinate 
transformations and gauge freedoms do not arise; this reduces the risk of direct 
computational mistakes, as well as the omission of special cases. It is an impor- 
tant development that the GHP formalism now has computer support [1], [2] 
and that attempts are being made to develop the GHP integration procedure 
algorithmically in the programmes [3] ; hopefully similar developments will occur 
soon for GIF. 

In [7] it was easy (in the generic case) to draw conclusions about Killing vector 
structure and invariant classification; this was because only intrinsic elements 
of the GIF were used. In this paper we have introduced a non-intrinsic element 
early in our calculations, and this makes it more difficult to draw conclusions 
about Killing vector structure and the invariant classification. We shall discuss 
in detail how to incorporate these consideration in general into GIF in a sub- 
sequent paper. At a more ambitious level, we are also applying this method to 
more general vacuum type N metrics. 
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